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ABSTRACT 

Coherent 65 mHz pulsations in the X-ray flux of the Small Magellantic Cloud (SMC) transient source 
RX J0052. 1-7319 have been detected by us in an analysis of ROSAT data. We report on the pulsations 
we detected in ROSAT HRI data and simultaneous detection of these pulses in hard X-rays using 
BATSE data. The BATSE data show an outburst of the source lasting 60 days. We report on optical 
observations of the candidate companion, and a new source position we determined from the HRI data, 
which is consistent with the candidate's location. From the measured fluxes and observed frequency 
derivatives we exclude the possiblity that the pulsar is in the foreground of the SMC, and show that an 
accretion disk is present during the outburst, which peaked near Eddington luminosity. 

Subject headings: accretion, accretion disks - binaries: general - pulsars: individual (RX J0052. 1-7319) 
- X-rays : stars 



1. INTRODUCTION 

RX J0052. 1-7319 is an X-ray source located in the the 
Small Magellanic Cloud, first detected with Einstein (IE 
0050.3-7335, Wang & Wu 1992). It was classified by Ka- 
habka and Pietsch (1996) as a transient X-ray binary can- 
didate based on the ROSAT PSPC data from October 
1991 and April 1992. 

As part of a systematic search of the ROSAT data for 
pulsed sources, we have discovered pulsations from RX 
J0052. 1-7319 at a frequency of 65.4 mHz (period = 15.3 
s) using ROSAT HRI observations from 1996 November- 
December. Our search had previously detected another 
SMC pulsar, J0117.6-7330 (Macomb et al. 1999). 

After a preliminary report of our discovery of the pulsar 
nature of RX J0052. 1-7319 (Lamb et al. 1999), Kahabka 
estimated a new source position based on HRI observa- 
tions of 1995 May (Kahabka 1999a), and confirmed the 
detection of pulsations using an HRI observations from 
October 1996 (Kahabka & Pietsch 1996). Israel et al. 
(1999) identified a Be star as the likely optical counterpart 
to the source. Udalski (1999) then reported on the long 
term optical variability of this candidate, based on Optical 
Gravtational Lensing Experiment (OGLE) monitoring. 

Here we report on the pulsations we detected in the 
ROSAT HRI data, detection of these pulses in hard X- 
rays using simultaneous BATSE data, the history of the 
outburst visible in the BATSE data, and optical observa- 
tions of the candidate companion and nearby stars. 

2. OBSERVATIONS 

'^ Astrophysics Programs, Universities Space Research Association 



We initially discovered pulsations from RX J0052. 1-7319 
in observations made with the High Resolution Imager 
(HRI) detector of the ROSAT (Triimper 1983) in 1996 
November and December. The HRI, which consisted of 
two cascaded microchannel plates (MCPs) with a crossed 
grid position readout system, was sensitive to X-rays in the 
~ 0.2 - 2 keV range. Thereafter we detected the pulsa- 
tions in observations from the Burst and Transient Source 
Experiment (BATSE) (Fishman 1989) on the CGRO us- 
ing data from the Large Area Detectors (LADs) which 
are Nal scintillation detectors sensitive to hard X-rays/soft 
gamma-rays in the 20 keV to 2 MeV range. 



2.1. ROSAT Observations 

Our discovery was made in the ROSAT HRI observa- 
tions of 1996 November 10.71 - December 9.16 as part of 
a systematic search of the ROSAT data for previously un- 
detected pulsars. These data were among 1365 data sets 
we selected from the catalog of 59911 ROSAT/HRI point 
source observations based on the potential for detection of 
significant pulsations. Each data set was processed using 
standard Ftools, with barycentered arrival times for the 
first 200 ks of each observation binned with 5 ms resolu- 
tion, and Fourier transformed. The power spectra were 
searched for significant pulsed signals unassociated with 
the ^5760 s spacecraft orbit or the 402 s period spacecraft 
orientation wobble. The RX J0052. 1-7319 HRI observa- 
tions were selected for further analysis because of a strong 
signal near 65 niHz. 
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Fig. 1. — Mean count rates in 1000 s bins for the 1996 Novem- 
ber 11 (MJD 50398) - December 9 ROSAT HRI observations of 
RX J0052. 1-7319 showing the time structure of the observation 
window, and the declining flux of the source. The upper panel 
is for observation rh600811n00 and the lower for rh600812a01. 
Between the panels is shown both the time interval for the initial 
FFT in which the pulsations were discovered, and the interval 
used for the analysis resulting figures 2 and 3. 

These HRI observations, rh600811n00 and rh600812a01, 
are of two SMC fields which contain RX J0052. 1-7319. The 
observations for these two fields are partially contempora- 
neous. Figure 1 shows the count rates in 1000 s bins for 
events within 60" of the source in these observations. The 
roughly 10-20% difference in counting rates between con- 
temporaneous portions of the two observations may be due 
principally to differences in source vignetting. For obser- 
vation rh600812a01, the source is 17' from the center of the 
HRI field and vignetting corrections are 10-20%; for ob- 
servation rhGOOSllnOO the source is 8'4 from the HRI cen- 
ter and vignetting is substantially less. From observation 
rh600811n00 it is clear that the source fiux is declining. 

As can be seen in Figure 1, the initial 200 ks (2.3 days) 
interval of the observations, which was used for the initial 
search fft, contained only a portion of either of the two 
observations. We chose for further analysis the data of 
Nov 11.21-16.83 (aU of rh600812a01 and the first portion 
of rh600811n00). 

From initial analyses of this dataset we knew pulse fre- 
quency was changing significantly, and that there was 
substantial power at higher harmonics of the pulse fre- 
quency. To accurately estimate the pulse frequency and 
frequency derivitive we maximized the Z| statistic (Buc- 
cheri et al. 1983), where Z|(/) = p{f) + p{2f) +p{3f), 
and the Rayleigh statistic p{f) is calculated as: 
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Fig. 2.— Rayleigh periodogram for the 1996 Nov 11.21-16.83 
ROSAT HRI observations of RXJ 0052.1-7319. This was imple- 
mented with an FFT of data binned with 0.5 s resolution. A 
fractional frequency derivative of v/v = 8.37 X 10"^" s~^ is ac- 
counted for in the analysis. Only values above 10 are plotted. 
The pulse frequency, near 65 niHz, and four harmonic overtones 
are clearly detected. The structure in the wings of each peak is 
due to the non-uniform spacing of the data. 

Here N is the number of photons detected in the inter- 
val, / is the analysis frequency, a — v jv^ with v the pulse 
frequency, ifc is a barycentric arrival time of photon fc, and 
T is an epoch within the data interval. 

For epoch MJD 50401.0 (1996 Nov. 14.0) we obtained v 
= 0.06545850(7) Hz and i> ^ 5.48(7) x 10"" Hz s^^. The 
maximum value of Z| was 524, which is highly improbable 
by chance.^ The Rayleigh statistic for this ratio of z>/j/ is 
shown in figure 2. In addition to the pulse fundamental, 
four pulse harmonic overtones are clearly detected. 

The HRI data from the 1996 November 11.21-16.83 in- 
terval epoch-folded with this ephemeris is shown in figure 
3. The pulse profile has a single asymmetric peak, with a 
narrow valley at minimum. The pulsed fraction ([mean- 
minimum] /mean) is 35 ± 3%. 

For the HRI data from the interval 1996 December 5-9 
we obtain a frequency i' = 0.0655330(3) Hz at epoch MJD 
50524.0, and i> = 2.3(2) x 10"" Hz s'^ Thus the de- 
crease in counting rate noted above was accompanied by a 
decrease in spin-up rate, as would be expected from the in- 
trinic correlation between mass accretion rate and angular 
momentum accretion rate in a disk fed accreting pulsar. 

For the 1996 November 11-16 and December 5-9 inter- 
vals we find source count rates of 0.92 and 0.46 counts s~^ 
respectively. These rates are corrected for vignetting 
(David et al. 1999) and the 15" radius circle used to select 
the data, using the HRI off-axis encircled energy profile 
(Boese 2000). 



^We estimate that the probability of obtaining a Z| of 524 or more due to Poisson noise is less than 10"^*, including the number of trials 
introduced by searching in frequency and frequency derivative. This calculation however neglects systematic signatures in the data, which 
would likely play a dominate role any false detection with a Z| this large. 
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Fig. 3.— Pulse profile of RXJ 0052.1-7319 from the ROSAT 
HRI observations of 1996 Nov 11.21-16.83. 

To estimate fluxes from these count rates we have used 
the spectrum estimated for the source by Kahabka & Pri- 
etsch (1996) based on a fit of ROSAT PSPC data. This 
consists of thermal brcmsstrahhing model with kT =17 
keV, a galactic absorption column of 3 x lO^^cm^^, and 
an absorption column within the SMC with density of 
Nh = 7.0 ± 4.0 X 10^-^cm~^ (95% confidence errors), with 
metallicities reduced by a factor of seven from solar sys- 
tem abundances. This model results in 0.1-2.0 keV unab- 
sorbcd fiux estimates of 5.9 ± 1.0 x 10~^^erg cm~^s~^ and 
of 3.0 ± 0.5 X lO^^^erg cm^^s^^ for the 1996 November 
11-16 and December 5-9 intervals respectively. 

Using the 1996 November data of observation 
rh600811n00, we have determined the position of the 
source to be R. A. = 0'^52™13!65, Dec. = -73°19'19'.'5, 
with a statistical uncertainty (1") which is negligible in 
comparison to systematic position errors, which result in 
an 68% confidence error circle of radius 7" (Kiirster 1993, 
ROSAT Users Handbook^). These systematic errors are 
due to uncertainty in the aspect of ROSAT. We note that 
this position differs by nearly 10" from the position quoted 
by Kahabka (1999b) and Kahabka (2000). However, since 
these latter determinations used either very low rate data 
(<0.01 c/s) or data in which the source was at the edge of 
the HRI field of view, we believe that the position given 
above may be the more accurate. 

2.2. BATSE Observations 



to 2000 May, except during SAA passages and occasional 
telemetry outages. 

The 20-50 keV channel of the BATSE LAD discrimi- 
nator rates (DISCLA channel 1) were analysed for pulsa- 
tions from RX J0052. 1-7319 using techniques discussed in 
Finger et al. (1999). Rates were combined from differ- 
ent detectors with coeflficients optimal for a source with a 
spectrum of the form dN/dE = A* exp{—E/kT)/E with 
kT — 20 keV. A large number of pulse profiles are obtained 
by fitting short segments of these combined rates with a 
model consisting of a quadratic spline background, plus 
a low order Fourier expansion pulse profile model. These 
profiles are then combined over multiple day intervals us- 
ing trial frequencies and frequency rates. The resulting 
combined profiles are evaluated with the Yn statistic (Fin- 
ger et al. 1999), which is simular to the Z^ statistic, but 
accounts for possible non-Poisson noise. 

Using data from the interval 1996 Nov 11.0-17.0, Y3 
was searched over a frequency range of width 10^* Hz 
centered on the ROSAT measurement, and a frequency 
rate range of —10^^^ Hz s"-'^ to lO""'^" Hz s~^ resulting 
in maximum value of ¥3 ~ 57.0 which we estimate"^ has 
a probability of being exceeded by chance of less than 
2 X 10~^. The resulting frequency and frequency rate esti- 
mates were v = 0.06545830(10) Hz at epoch MJD 50401.0, 
and i> = 5.62(12) x 10^^^ Hz s~^, in good agreement with 
the ROSAT results. 

The BATSE pulse profile for this interval using the 
ROSAT HRI ephemeris is shown in figure 4. The pro- 
file is shown relative to the mean flux level, which can- 
not be determined from the data due to the weakness 
of the source and the high background level. The solid 
curve is the proflle corresponding to six Fourier coefficients 
which have been estimated with the same fitting technique 
used with the frequency search. Error bars for the value 
of the curve at 13 approximately independent points are 
shown. Six harmonics were chosen because this approxi- 
mately matches the 1.024 s resolution of the data. 

Pulsations were then detected in nine additions six-day 
intervals, extending the total period the source was de- 
tected with BATSE to 60 days between 1996 September 18 
and November 17. Beginning with the Nov 11-17 interval, 
the total period of detection was progressively extended 
outward by frequency and frequency rate searches in the 
adjacent six day intervals, using 5 x 10~^ Hz frequency 
ranges centered on the extrapolation from the neighbor- 
ing frequency and frequency rate estimate, and the fre- 
quency rate range given above. Searching in this manner 
at the boundaries of the known frequency history reduces 
the required frequency search range, providing the best 
sensitivity. Figure 5 shows the resulting measurements of 
frequency, frequency rate, and r. m. s. pulsed flux in the 
20-50 keV energy band. The figure also shows for compar- 
ison the frequency and frequency rate measurements from 
the ROSAT HRI data. 



The Burst and Transient Source Experiment (BATSE) 
was an all sky monitor which flew onboard the Comp- 
ton Gamma-Ray Observatory (CGRO). The data used in 
the analysis presented here were from the Large Area De- 
tectors (LADs) which were Nal(Tl) scintillation counters 
1.27 cm thick, with 2025 cm^ area, that were located at 
the eight corner of the CGRO spacecraft. Count rate data 
from the LADs was continuously available from 1991 April 

•^ http://heasearc.gsfc. nasa.gov/docs/rosat/ruh/handbook/nodc34. html 

^This estimate approximates the distribution of Y3 with the chi-squared distribution with 6 degrees of freedom, which is valid in this case 
because of the large number of profiles being combined, and accounts for 1100 independent trials. 
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reasonable counterpart is object A, the star identified by 
Isreal et al. (1999). Object A is well contained with the 
error circle for RX J0052. 1-7319 presented here, and thus 
the identification of the X-ray source as a member of a Be 
system seems reasonable. 

The absolute V magnitude of object A may be de- 
termined from the value quoted by Udalski (1999) of 
m^ = 14.67. The distance modulus to the SMC determined 
by Westerlund (1997) is (m - M)o = 18.9. In addition, 
the average extinction is E(B-V) — 0.07-0.09 (Schwering & 
Israel 1991), though there are regions in the SMC where it 
can rise as high as 0.25. Combining these parameters leads 
to an estimate of M^= -4.47±0.02. This magnitude corre- 
sponds to a star in the range BlIII - BOV, very similar to 
the value obtained for Be/X-ray binary counterparts (e.g. 
that of RX J0117.6-7330 quoted in Coe et al. (1998)). 
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Fig. 4.— Pulse profile of RXJ 0052.1-7319 from the BATSE 
observations of 1996 Nov 11.0-17.0. This profile is based on an 
estimated six harmonic Fourier model. Error bars are given for 
approximately independent points. This profile is folded with 
the same ephemeris at that in figure 3. 

Efforts were made to detect the source from 1996 Au- 
gust 1 to September 18, and from 1996 November 17 to 
1997 January 8, but these failed. These attempts used 
searches with a frequency range of 3 x 10^'^Hz centered 
on the nearest detection, and a frequency rate range of 
-10"" Hz s-i to 10-1° Hz s-^ Prior to 1996 Septem- 
ber 18, pulsations from the 130.4 mHz pulsar 4U 1626-67 
interfered with first harmonic (i.e. 2 x /) contribution to 
Ya, in a narrow frequency band near 65.2 niHz. For this 
narrow band, this harmonic was left out of the search. We 
conclude that in these intervals prior to and following our 
1996 September 18 to November 17 detections that the 
source had pulsed flux below our sensitivity level. We es- 
timate an upper limit of 8 mCrab for the 20-50 keV r. m. s. 
pulsed flux for these intervals. 
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2.3. Optical Observations 

In figure 6 we show a V band image taken from South 
African Astronomical Observatory (SAAO) 1.0m telescope 
on 20 Jan 1999. Marked on the figure is the X-ray error 
circle of Kahabka (2000) which has a 6" radius, and the 
error circle from this work, with a 7" radius. The object 
marked A is the source identified as the counterpart by Is- 
rael et al. (1999) based upon a red spectrum. The objects 
labeled A, B, C, and D are also identified in Figure 7. 

If the optical counterpart of RX J0052. 1-7319 is a Be 
star, then it should show excess H^ emission. In figure 7 
we show Ha versus R band counts for 25 stars in the vicin- 
ity of Kahabka's X-ray error circle. The line is a linear best 
fit through all the data points shown. Stars that lie above 
the line are ones that show an Ha excess - i.e. objects A, 
B and C as identified on figure 6. Object D, the one inside 
the Kahabka (2000) error circle, shows no excess. Objects 
B and C, though exhibiting an Ha excess, are too far away 
from the X-ray position to be the counterpart. If, as we 
might expect, the optical companion is a Be star, the only 
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Fig. 5. — Pulse timing measurement for RX J0052. 1-7319 
from BATSE and ROSAT, and pulsed flux measurements from 
BATSE. The top panel shovifs the pulse frequency, with solid cir- 
cles for BATSE measurements, and open diamonds for ROSAT 
HRI measurements. The lines through each point span the mea- 
surement interval, and have the slope of the measured frequency 
rate. The error bars are too small to be visible on the plot. 
The middle panel shows the pulse frequency rate, with solid cir- 
cles for BATSE measurements, and open diamonds for ROSAT 
HRI measurements. The lines through each point span the mea- 
surement interval. The bottom panel shows the BATSE r. m. s. 
pulsed flux (20-50 keV). For reference 1 mCrab (20-50 keV) = 
1.0 X 10~^^ erg cm^ s~^, which corresponds to a luminosity of 
4.3 X lO'^^erg s~^ in the SMC. The plot covers the date range 
1996 September 14 to December 23. For flux upper limits before 
and after the BATSE measurements presented see the text. 




Fig. 6. — V band image taken from South African Astronom- 
ical Observatory (SAAO) 1.0m telescope on 20 January 1999. 
Marked on this figure is the X-ray error circle reported in Ka- 
habka (2000) near Star D, and the error circle from this work 
located near Star A. 



3. DISCUSSION 

The ROSAT HRI observations we have presented show 
that RX J0052. 1-7319 is a pulsar with a 65 niHz rota- 
tion frequency. Variations in the flux by more than two 
orders of magnitude between ROSAT observations at dif- 
ferent epochs demonstrates that the source is a transient 
(Kahabka & Pietsch 1996,Lamb et aL 1999). The BATSE 
observations we have presented show that the November - 
December 1996 ROSAT HRI observations occurred on the 
tail of a outburst that began at least two months earlier, 
lasting more than 80 days. 

The companions of transient pulsars with known spec- 
tral type are generally Be (or Oe) stars. There are now 
known a few transient pulsars in low- mass X-ray binaries, 
but these all have high spin frequencies {ly > 1 Hz). Be 
stars show both emission in the Balmer lines, and an ex- 
cess in the IR, due to a disk of material shed from the 
equator of the rapidly rotating star. The accretion of this 
material, which is thought to form a quasi-keplerian disk, 
fuels the X-ray outbursts. 

Two types of outbursting behavior are observed in 
Be/pulsar binaries: type I ("normal" outbursts), with a 
series of lower luminosity (L < 10^^ erg s~^ outbursts 
occurring once per orbit, generally near periastron; and 
type H ("giant" outbursts), single high luminosity {L 



2^q38 gj.g g-i sometimes lasting several orbits (Stella, White 
& Rosner 1986). These outbursts are generally accompa- 
nied by rapid spin-up of the neutron star, indicating disk 
accretion (Finger, Wilson & Harmon 1996, Bildsten et al. 
1997). 

The BATSE observations presented here are consistent 
with a giant outburst of a Be/pulsar transient in the 
SMC. From the observed spin- up we can show that a disk 
is present, and the the source is outside our galaxy. In 
disk accretion, the angular momentum accreted per unit 

accreted mass is / = {GMr„i) '^, were M is the neu- 
tron star mass, and r^ the radius of the magnetosphere. 
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Fig. 7. — Ha versus R band plot of 25 stars in the vicinity of 
the X-ray error circle. Stars that lie above the line are ones that 
show an H alpha excess. Letters refer to stars identified in the 
previous figure. 

The specific angular momentum for wind accretion is gen- 
erally much less than this, since if it reaches this value, a 
disk forms. The magnetospheric radius is limited by the 

1/ 2/ 

corotation radius Tc — (GM) /3(27ri^)^ /s — 1.0x10^ cm, 
where we have used M ~ IAMq. The specific angular 
momentum / is proportional to the ratio ;//L where L is 
the pulsar's luminosity. Since the rate at which angular 
momentum is accreted is Irh = 2'kIv where I is neutron 
star moment of inertia, and L — GMm/ R, where R is the 
neutron star radius, we have the upper limit 



< (27r)"3(GM)" 
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= 3.8 X 10"'^-^ Hz erg" 



where we have used a moment of inertia / 
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"gm cm^ 
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and a radius R - 
vation, for which 
L > 1.5 X 10^8 



10 cm. For the Nov 11-17 obser- 



= 5.6 X 10^^ Hz s , this implies 
'^ and for peak frequency rate ob- 
Using the 0.1-2.0 
keV flux F — 5.9 x 10^^^ erg cm^^ s^^ observed in the 
1996 Nov 11-17 HRI observation, we find the lower limit 
on the source distance of 



erg s 
served we have L > 2.0 x 10^® erg s^^. 
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where e is the fraction of the luminosity in the 0.1-2.0 keV 
energy band. This eliminates the possibility that the RX 
J0052. 1-7319 is in the foreground of the SMC. 

The limit in equation 2 is reached only with disk accre- 
tion. Using the measured frequency rate, fiuxes, and the 
distance of 60 kpc to the SMC we find a ratio of i'/L com- 
parable to this limit, suggesting that an accretion disk is 
present during the outburst. 

The pulse periods of Be binary pulsars range from 69 ms 
to 1400 s, with determined orbital periods ranging from 17 
to 250 days. Corbet (1986) showed that the orbital period 
is correlated with the pulse period. From the observed 
distribution of spin and orbital periods we would expect 
an orbital period in the range of 25 - 100 days. In giant 
outbursts of Be/X-ray pulsars we expect a strong correla- 
tion of flux and frequency rate (see e.g. Finger, Wilson & 
Harmon 1996). We note however that the history of the 
frequency rate in fig 5 is dissimilar in profile to that of the 



pulsed flux. This could be due to the doppler signature of 
a binary orbit. An alternate explanation would be changes 
in the spectra or pulse fraction. However, giant outbursts 
typically have simple flux and intrinsic spin-up rate pro- 
files, with a steady rise to peak, and a somewhat slower 
faU (see e.g. Parmar et al. 1989, Whitlock 1989, Finger, 
Wilson & Harmon 1996, but also, Negueruela et al. 1997). 
The plausible identification of RX J0052. 1-7310 as Be 
binary system in the SMC accentuates further the rather 
dramatic difference between the SMC and our Galaxy with 
regard to the population of high mass X-ray binaries. This 
fact has already been noted by several authors (Schmidtke 
et ak, 1999, Jokogawa et al. 2000). A recent compila- 
tion of the known X-ray pulsars'* gives within the SMC 
one supergiant system, three known Be systems, and 11 
transients with uncertain companion class (likely to be Be 
systems), making 15 high mass X-ray pulsar binaries. For 
the Galaxy the corresponding number is 40. Therefore, 
using a mass ratio of the SMC to the Galaxy of 1/100, 
this suggests that high mass X-ray pulsar systems in the 
SMC are overabundant by roughly a factor of 30 relative 
to the Galaxy. This analysis ignores the important effects 
of obscuration within the Galaxy, the relative frequency of 
observations, and the low luminosity sensitivites obtained 
for the Magellantic clouds; nevertheless the apparent dis- 



parity is remarkable. 

Since high mass X-ray binaries have lifetimes which are 
a very small fraction (~ 10""^) of the age of the Galaxy, 
the dramatic difference between the SMC and the Galaxy 
points to a rather recent outburst of star-formation in the 
SMC within the last ~ 10^ years. Further support of such 
an epoch of star formation comes from the radio observa- 
tions of HI by Stavely-Smith et al. (1997) and Putman et 
al. (1998), which show a strong bridge of material between 
the Magellanic Clouds and between them and our own 
galaxy. Furthermore, Stavely-Smith et al. have demon- 
strated the existence of a large number of supershells (cre- 
ated by multiple supernovae) of a similar age (~5 Myr), 
strongly suggesting enhanced starbirth has taken place as 
a result of tidal interactions between these component sys- 
tems. Consequently it seems very likely that the previous 
closest approach of the SMC to the LMC ~ 10^ years ago 
may have triggered the birth of many new massive stars 
which have given rise to the current population of HMXBs. 
In fact, other authors (e.g. Popov et al. 1998) claim that 
the presence of large numbers of HMXBs may be the best 
indication of starburst activity in a system. 

M. H. F. acknowledges support through NASA grant 
NAG5-4238. 
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